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Abstract 

We searched for quantitative trait loci (QTL) associated with the palm oil fatty acid composition of mature fruits of the oil 
palm £. guineensis Jacq. in comparison with its wild relative £ oleifera (H.B.K) Cortes. The oil palm cross LM2T x DAIOD 
between two heterozygous parents was considered in our experiment as an intraspecific representative off. guineensis. Its 
QTLs were compared to QTLs published for the same traits in an interspecific Elaeis pseudo-backcross used as an indirect 
representative of £ oleifera. Few correlations were found in E. guineensis between pulp fatty acid proportions and yield 
traits, allowing for the rather independent selection of both types of traits. Sixteen QTLs affecting palm oil fatty acid 
proportions and iodine value were identified in oil palm. The phenotypic variation explained by the detected QTLs was low 
to medium in £ guineensis, ranging between 10% and 36%. The explained cumulative variation was 29% for palmitic acid 
CI 6:0 (one QTL), 68% for stearic acid CI 8:0 (two QTLs), 50% for oleic acid CI 8:1 (three QTLs), 25% for linoleic acid CI 8:2 (one 
QTL), and 40% (two QTLs) for the iodine value. Good marker co-linearity was observed between the intraspecific and 
interspecific Simple Sequence Repeat (SSR) linkage maps. Specific QTL regions for several traits were found in each mapping 
population. Our comparative QTL results in both £ guineensis and interspecific materials strongly suggest that, apart from 
two common QTL zones, there are two specific QTL regions with major effects, which might be one in £ guineensis, the 
other in £ oleifera, which are independent of each other and harbor QTLs for several traits, indicating either pleiotropic 
effects or linkage. Using QTL maps connected by highly transferable SSR markers, our study established a good basis to 
decipher in the future such hypothesis at the Elaeis genus level. 
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Introduction 

Indigenous to Africa, the oil palm {Elaeis guineensis Jacq.) is a 
perennial, monocotyledonous, monoecious, cross-pollinating spe- 
cies belonging to the Arecaceae family. The only other species in the 
genus Elaeis is the American oil palm, Elaeis oleifera (H.B.K) Cortes, 
indigenous to the Amazon region in South America [1,2]. Both 
species have 16 chromosome pairs (2n = 32) [3], and they can 
easily hybridize with each other [4] . 

Beginning the second year after planting and continuing 
throughout its life, the cultivated oil palm produces unisexual 



male or female inflorescences in successive cycles, emerging at the 
axU of each leaf Female inflorescences grow in bunches that hold 
between 200 and 4 000 fruits. The oil palm fruit is a drupe. It 
comprises a pulp (mesocarp), an endocarp, called the shell; and a 
kernel. Three fruit types exist, depending on the presence or 
absence of the shell, which is governed by a major gene called Sh 
[5] . The dura type, homozygous Sh^/Sh^, produces large fruits with 
a thick shell and a pulp that is fairly abundant by weight (35-70%). 
The pisifera type, homozygous Sh~ /Sh~ , is generally female-sterile, 
and its few fruits are relatively small with a relatively large pulp 
(90%). The tenera type, heterozygous genotype SH^/Sh , thin 
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shelled with an abundant pulp, produces the most palm oil and is 
therefore the fruit type of aU commercial oU palm varieties. Only 
the thick-shelled dura type exists in the E. oleifera species. 

The oil palm produces two distinct vegetal oils in its fruits: crude 
(red) palm oil, the fresh oil obtained from the mesocarp, and kernel 
oil, the oil obtained from the kernel seed of the palm fruit. These 
oils are used and marketed separately according to their own 
supply and demand conditions [6] . The oil productivity of the oil 
palm is ten-fold higher than that of the soybean, with a yield 
reaching 7.5 to 9 tons of oU/year/ha for the best Cirad 
commercial varieties in favorable agro-climatic conditions (www. 
paJmeUt.com). Palm oil constitutes the highest annual world 
production of vegetal oil, at 53.3 Mt, followed by soybean oU 
(43.4 Mt) and rapeseed oil (23.5 Mt) (http://www.fas.usda.gov/ 
oilseeds/Current/). The wild species E. okifera is not exploited at a 
commercial level due to its low yield [7] . 

The consumption of vegetal oils falls into two major applica- 
tions: the food industry (with over 80% of the market) and the 
chemical industry [8]. Furthermore, in the 1980s, with the start of 
the full exploration of vegetable oils for biodiesel, principally due 
to interest in renewable energ)^ sources, oil crops such as oil palm, 
soya, and oilseed rape began to be used as primary sources for 
these fuels, and higher proportions of unsaturated fatty acids are 
better for this purpose [9,10]. 

Palm oil contains approximately 50% saturated fatty acids, with 
44% pahnitic acid (C16:0), 5% stearic acid (C18:0), and trace 
amounts of myristic acid (C14:0). The unsaturated fatty acids are 
approximately 40% oleic acid (CI 8:1) and 10% polyunsaturated 
linoleic acid (C18:2) and linolenic acid (C18:3) [11-13]. 

Currendy, the need for high-oleic crops is increasing as the food 
market and the agro-industry demand oils that are more resistant 
to oxidation. Because of its saturated/unsaturated fat ratio close to 
1 and its concentrations of carotenoids, tocopherols, and 
tocotrienols, palm oil is considered an oxidatively stable oil 
[6,14] that has no equivalent among other vegetal oUs. Another 
goal linked to the high oleic acid content is to improve the frying 
properties of paJm oil for the food market, although the quality of 
palm oil in deep-fat frying methods has already been established 
[15]. Consequently, the improvement of palm oil quality for a 
higher degree of fatty acid unsaturation is one of the goals of oil 
palm breeders. 

Several studies have investigated the phenotypic characteriza- 
tion and genetic determinants of the palm oil fatty acid 
composition [12,16-18]. These studies showed that there is 
significant variability in the fatty acid composition among E. 
guineensis populations. In addition, the composition of palm oil is 
substantially different between E. guineensis and E. okifera, the latter 
characterized by its higher content of unsaturated fatty acids. 
Among the E. oleifera palms, the unsaturated fatty acid content 
ranges from 47% to 69% for C18:l, 2% to 19% for C18:2, and 
0.1% to 1.2% for C18:3 [19-21]. Consequently, the iodine value 
(IV), a multi-parameter measure of the global degree of 
unsaturation of the fatty acids in vegetal oU, is higher in E. okifera 
than in E. guineensis due to the major oleic acid C18:l and the 
linoleic acid C 1 8:2. The IV for E. oleifera is between 70% and 87% 
[13,19,20,22], whereas the value for E. guineensis varieties is 
between 53% and 60% [12,22-24]. Studies in interspecific hybrids 
showed that the palm oil composition profile is intermediate 
between those of the two parental Elaeis species for most fatty 
acids, which is an advantage for introgressing these traits from E. 
okifera into varieties of E. guineensis [4,7,25] . As far as we know, only 
two QTL studies have been performed for palm oil fatty acid 
composition, one in oil palm [60] , the other one in an interspecific 
pseudo-hybrid [26], and no published study has compared QTLs 



of the palm oil fatty acid composition between the two Elaeis 
species or between intraspecific and interspecific materials except 
for Montoya et al. [27]. 

According to the review by Murphy [28] and references therein, 
the oil crop varieties that have acquired significantly increased 
unsaturated fatty acid proportions by conventional breeding 
methods are the rapeseed/canola, soybean, sunflower and 
safilower (with 75% oleic acid and 1% to 3% linolenic acid), 
and the olive (with 85°/() oleic and l""!) linolenic acid). Additionally, 
transgenic approaches have produced varieties that have signifi- 
cantly increased unsaturated fatty acid content: rapeseed/canola, 
with 89% oleic acid; Indian mustard, with 73% oleic acid; 
soybean, with 75%o oleic acid; and cotton seed, with 78% oleic 
acid. One breeding strategy to achieve a highly unsaturated pabn 
oil is to take advantage of the unsaturated profile found in E. 
olefera using interspecific hybridization or backcross strategies 
[4,7,29]. 

Few authors have reported data describing correlations, at the 
level of individual pakn trees between fatty acid proportions in the 
Elaeis genus. Even fewer data have been published regarding the 

correlation of fatty acid proportions with vegetative and produc- 
tion traits in support of breeding strategies to modify the palm oil 
composition. 

By contrast, in this respect, there is evidence of interference 
between agronomic traits and oil quality, as reported for the 
soybean. Severed studies have shown that the presence of major 
and minor genes that reduce the palmitate (C 16:0) content in seed 
oil also reduces seed yield or plant height due to pleiotropic effects 
or linkage with unfavorable yield genes [30,31]. However, the 
increased oleic acid in soybeans with high-yielding genetic 
backgrounds does not aflect yield or other agronomic traits 
[32,33]. 

In the olive tree [Olea europaea L.), correlations between fruit 
characteristics, oil yield components, and fatty acid composition in 
j)r()gc'nies from different crosses show significant positive relation- 
ships between oil content and oleic acid percentage, which are 
negatively correlated with the palmitic, palmitoleic, and linoleic 
acid contents [34]. 

The principal aim of the present study was to identify the 
genetic determinants of the palm oil fatty acid composition in the 
cultivated oil palm in comparison with its wild relative E. okifera. 
As a preliminary study, the palm oil composition was character- 
ized in an intraspecific E. guineensis cross, and the relationships 
between its components were studied to better understand the 
genetic determinants of palm oil fatty acid proportions. In 
addition, relationships with other production traits were studied 
to determine whether selection for a modified palm oil compo- 
sition could significantiy influence the pakn oil yield traits. Based 
on the phenotypic variability existing between E. guineensis and E. 
olefera in the profile of fatty acids of their respective oils, we 
hypothesized that the genetic architecture of this trait would difler 
between these E. guineensis intra-specific and interspecific mapping 
populations, in terms of the number and position of QTLs and in 
the phenotypic variance they explain. To test this hypothesis, two 
Elaeis mapping populations were exploited: an E. guineensis 
intraspecific cross and an interspecific pseudo-backcross. A second 
hypothesis was that the two Elaeis genomes would have good locus 
coUinearity because they are not distant relatives, as they intercross 
easily. In these crosses, dense microsateUite linkage maps, with a 
high number of common simple sequence repeat (SSR) marker 
loci, allowed us to identify and compare QTLs. The results of the 
QTL detection undertaken in the E. guineensis cross were 
interpreted in comparison with QTLs for the same traits in the 
interspecific pseudo-backcross. 
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Previous QTL detections for palm oil fatty acid composition in 
the Elaeis genus [26,60] where unfortunately based on incomplete 
linkage maps mostly using AFLP, RFLP, or SSCP non-transfer- 
able marker loci. Such situations did not allow to relate and to 
compare QTL results, as well as with ours and those of Montoya et 
al. [27] which are in the contrary based on a saturated linkage 
map, made of SSR marker loci higly polymorphic and transferable 
within the EMeis species. Therefore, we chose to use same linkage 
map quality and related E. guineensis map parents in our 
comparative experiment. 

Materials and Methods 

Vegetal material 

The E. guineensis cross LM2T x DAIOD, previously described by 
Billotte et al. [35], was used to perform a QTL analysis of fatty acid 
proportions in palm oil from mature fruits. The progeny consisted 
of 116 fuU-sibs derived from the cross between two heterozygous 
E. guineensis parents from the CNRA oil palm breeding program 
(La Me, Cote d'lvoire): LM2T, a tenera palm belonging to the La 
Me African population, and DAIOD, a dura palm selected from the 
Deli population. In this cross, traceable DNA marker alleles (and 
detectable QTL marker alleles) are those segregating from 
heterozygous loci of one, the other or both E. guineensis map 
parents (Figure SI A). 

We compared the intraspecrfic QTLs identified in LM2T x 
DAIOD to those shown to be responsible for the same traits in an 
Elaeis interspecific pseudo-backcross, SA569, described by Mon- 
toya et al. [27]. The SA569 mapping population consisted of 134 
fuU-sibs derived from an interspecific palm, SA65T [E. oleifera 
SA49D x E. guineensis LM2466P), and an E. guineensis genitor, 
P03228D. The female grandparent (SA49D) was a wild E. oleifera 
palm [dura) from the Coari region (Brazilian Amazon), and the 
male grandparent (LM2466P) was a pisifera oil palm obtained by 
selfing the LM2T oil palm genitor (parent of the intraspecific E. 
guineemis cross mentioned above). The E. guineensis male parent 
P03228D was a dura oil palm derived by selfing the DAI 15D oil 
palm genitor of the Deli population. Sharing same Deli ancestors, 
DA115D (subsequently its self P03228D) is genetically very close 
to the parent DAIOD of the E. guineemis cross mentioned above. 
The E. guineensis grand-parent (LM2466P) and parent (P03228D) 
are highly homozogous, over 75% [27]. By genetic construction of 
the pseudo-backcross, most traceable DNA marker loci are 
segregating only from heterozygous loci in the interspecific parent 
SA65T, each of them holding one E. guineensis allele (from LM2T 
or DA115D closely related to DAOD) and one E. oleifera allele 
(Figure SIB). Most QTL would be polymorphic and detectable on 
the interspecific parent SA65T only, while most of these QTL wiU 
be monomorphic and not identifiable on the E. guineemis parent 
P03228D, but detectable in the LM2T x DAIOD cross when 
existing and heterozygous on the E. guineensis LM2T or DAIOD 
parent. 

Elaeis linkage map system for intra- and interspecific QTL 
comparison 

1) Intraspecific multi-parent SSR consensus map in E. 

guineensis. This map, published by Billotte et al. [36] and 
henceforth referred to as Eg_Map, was constructed using SSRs 
genotyped in a 2 x 2 complete factorial mating experiment of four 
unrelated parents belonging to the La Me population from Africa 
(tenera LM2T), the Yangambi population from Africa (tenera 
LM718T), and the DeH population (dura DAIOD, dura LM269D). 
Each family was a single cross between two heterozygous parents, 
including one tenera parent from Africa and one dura parent of 



the DeU population: LM2T x DAIOD (mentioned above), LM2T 
X LM269D, LM718T x DAIOD, or LM7I8T x LM269D. 

These four crosses used by Billotte et al. [36] to establish the 

consensus Eg_Map were part of a larger genetic trial begun in 
1986 by the company SOGFINDO (Medan, Indonesia), whose 
experimental design was a randomized complete block design 
(RCBD) with five replications of 15 palms. The agro-climatic 
conditions were highly favorable for oil palm growth. The 116 
LM2T X DAIOD fuU-sibs described above were part of that 
experiment, and all other crosses were represented by 6 1 palms (a 
total of 299 palms for Eg_Map). The consensus Eg_Map, 
including 16 linkage groups, had 253 loci (251 SSRs, the Sh locus, 
which controls shell thickness, and its marker E-Agg/M-CAA132) 
and measured 1 731 cM (Haldane distance), with an average 
marker density of 7 cM [36]. 

2) Elaeis interspecific linkage map. The second genetic 
map used in this study was a dense SSR-based linkage map made 
from the interspecific pseudo-backcross SA569 described above, as 
published by Montoya et al. [27]. In this pseudo-backcross, 
molecular marker alleles were traced that segregated from both E. 
oleifera and E. guineensis grandparents and then from the interspe- 
cific and E. guineensis parents. The SSR map of SA569 had 362 loci 
(347 SSRs, 14 SNPs, and tiie Sh locus) and measured 1 485 cM 
(Hcddane distance), with an average marker density of 4 cM. In 
total, 156 marker loci (155 SSRs + Sh) were shared in common by 
E. oleifera and E. guineemis, and there was good marker co-linearity 
with Eg_Map, enabling the comparison of intra- and interspecific 
QTLs for fatty acid composition. 

Measurements of the palm oil fatty acid composition in 
LM2T X DAIOD 

The LM2T x DAIOD progenies, as well as five LM2T self 
palms and five DAIOD self palms, were analyzed for fatty acid 
composition and for the iodine value (IV) of the palm oil in mature 
fruits. The trait mean value from the LM2T and DAIOD self 
palms gave an estimate of the LM2T and DAIOD parent trait 
values. Of the 116 progenies planted initially for this cross, only 88 
were available that had produced bunches and allowed for 
estimates of the posterior palm oil fatty acid composition and 
iodine value due to the death of some palms before our phenotypic 
characterization in 2012. 

Measurements of the pakn oil fatty acid composition in the 
LM2T X DAIOD cross were performed as per Montoya et al. [27]. 
Quantitative phenotypic traits in this study considered the nine 
main fatty acids: myristic acid (CI 4:0), palmitic acid (C16:0), 
palmitoleic acid (CI 6:1), stearic acid (C18:0), oleic acid (CI 8:1), 
linoleic acid (C18:2), linolenic acid (C18:3), arachidic acid (C20:0), 
and gadoleic acid (C20:l). The iodine value (fV) was determined 
by the Wijs method described in the ISO 3961:2009 standard. For 
the IV, mean values were estimated for all 116 palms (based either 
on four measurements recorded in 2002-2004 and 2012 for the 88 
surviving palms or on two previous measurements in 2002-2004 
for palms that had died since 2004). 

Regarding the interspecific pseudo-backcross SA569, a total of 
115 progeny, as well as the E. olefera grandparent SA49D, the 
interspecific parent SA65T, and the E. guineemis parent P03228D, 
had been previously analyzed by Montoya et al. [2 7] for fatty acid 
composition and the iodine value of the palm oil in mature fruits. 
Available individual phenotypic data were used here for a 
comparative study with the LM2T x DAIOD cross. 
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Other phenotypic traits and QTL data available for LM2T 

X DA10D and SA569 

The fruit variety and the individual phenotypic data of 26 
vegetative or yield quantitative traits were available from Billotte 
et al. [36] for the four crosses of Eg_Map, including LM2T x 
DAIOD (Table SI). The main quantitative variables of production 
that we considered were the average bunch number/ palm/ year at 
3-5 years (Bn3_5), average bunch weight at 3-5 years (kg) 
(Bwt3_5), average bunch number/palm/year at 6-9 years 
(Bn6_9), average bunch weight at 6-9 years (kg) (Bwt6_9), average 
number of spikelets per bunch (Spikelets), average number of fruits 
per bunch (Fn), fruit to bunch ratio (%FB), average fruit weight (g) 
(Fwt), pulp to fruit ratio ('/oPF), palm oil to pulp ratio (%POP), 
kernel to fruit ratio (%KF), and iodine value (IV). Except the fruit 
type, no vegetative or production traits were available for the 
pseudo-backcross SA569. 

The QTL information pubhshed by Billotte et al. [36] for the 26 
vegetative and production traits and for the iodine value (IV) 
associated with the palm oil fatty acid profile in E. gumeensis (QTL 
position, confidence interval) were available for Eg_Map and 
therefore could be compared to the QTL information for pahn oil 
fatty acid composition identified for LM2T x DAIOD using 
Eg_Map. Based on the SA569 map, 19 QTLs for palm oil fatty 
acid composition and 14 Elaeis intra-gene SNP markers for five 
gene functions associated with oleic acid CI 8:1 in palm oil were 
available from Montoya et al. [27]. 

Statistical analysis of phenotypic data 

The iodine value was calculated from the average of all 
measurements made per pakn, as described above. Statistical 
analyses were performed on phenotypic data for fatty acid 
composition in the LM2T x DAIOD cross. The Gauss distribution 
of the quantitative data of the cross was checked by the Shapiro- 
Wilk normality test at the a threshold of 5 % . The relationships 
between phenotypic traits published by Billotte et al. [36] and the 
fatty acid proportions estimated in the LM2T x DAIOD cross were 
estimated by calculating Pearson's correlation coefficients at the 
individual palm level (Table S2). 

In a second step, the vegetative and production phenotypic data 
of LM2T x DAIOD were standardized as the m(;an and variance 
for both dura and tenera varieties to eliminate the Sh major gene 
effect on traits. The fatty acid composition data were not 
standardized because very few, neghgible correlations existed with 
the fruit type. All these data were subjected to an overall Pearson's 
correlation analysis to determine, at the individual palm level, the 
relationships between the fatty acid proportions and the vegetative 
and production traits. 

Principal component analysis (PCA) provides a synthetic vision 
of the relationships between studied variables based on estimated 
correlations. Every principal component (PC), a linear combina- 
tion of variables, is mathematically orthogonal to every other PC. 
Different PCs illustrate the degree of independence between their 
respective contributing variables. PCA was carried out using 
XLSTAT statistical software (Addinsoft, NY, USA) to visualize the 
associations among yield and palm oil composition traits, at the E. 
guineensis intraspecific level, using the Pearson's correlation matrix 
of LM2T X DAIOD. 

Associations among palm oil fatty acid composition traits were 
analyzed at both the intraspecific and interspecific levels by 
performing a PCA of five major fatty acids (CI 4:0 C16:0, CI 8:0, 
CI8:1, and C18:2) in LM2T x DAIOD and in SA569. 



QTL analysis of LM2T x DAIOD 

QTL analyses of palm oil fatty acid composition were 
performed with MapQTL5 [37] using the E. guineensis SSR 
linkage map (Eg_Map) and all available genot)rpic and phenotypic 
data for LM2T x DAIOD. No QTL analysis of other vegetative or 
production traits was performed here, as related QTLs are already 
available in Billotte et al. [36] . 

Three methods were used for QTL detection as per Montoya 
et al. [27]: i) A non- parametric Kruskal-Wallis (K-\\^ test was 
performed to identify significant marker-trait associations at p< 
0.005. ii) The interval mapping (IM) method was performed with a 
mapping step size of 1 cM and a maximum of five neighboring 
markers. To declare the presence of a QTL, the threshold LOD 
values were estimated at the genome-wide (GW) global risk a of 
5% and 1% by the re-sampling method and permutation of the 
trait data (1000 iterations), iii) The muItipIe-QTL model (MQM) 
method was carried out in conjunction with the automatic 
selection of cofactors, using the threshold LOD values described 
above. The threshold value of p<0.005 in the non-parametric K- 
W method was adopted on a empiric basis, as it was equivalent in 
QTL results compared to a GW global risk a of 5% using the 
different IM and MQM methods. A GW global risk ot of 1 % was 
also used in the IM or MQM method to explore the robustness of 
detected QTL. 

The confidence interval of each significant QTL by IM or 
MQM was determined by the LOD -1 method. For the 
phenotypic values of fatty acids (C16:l, C18:3, C20:0 and 
C20:l) in trace or in small amounts that did not follow a mixture 
normal distribution (Figure S2), the Kruskal-Wallis rank sum test 
was only considered as applicable for data with distributions far 
from mixture normal distribution [38]. A limited population size 
for identifying QTLs affects the accuracy of determining QTL 
locations and estimating QTL effects and, consequently, overes- 
timates the phenotypic variances associated with QTLs [39-41]. 
To correct at least the small part of the overestimation due to 
sampling error, we applied the correction described in Montoya et 
al. [27], as proposed by Luo et al. [42] and Xu [43]. They suggest 
to multiply the explained variance by 1 — l/(2*Ln(10)*LOD). 

Therefore, considering that the variance explained by an 
identified QTL is, as estimated under MapQTL, 



% variance explained = 100 ^cr^/cr^j 

with (j^ corresponding to the genetic variance due to additive 
effect and Up^ corresponding to the phenotypic variance, 

the corrected variance explained by this QTL was re-estimated 
like follows: 

% Corrected variance explained 
= (c7aVc7p2)(l-l/2Ln(10)xLOD) 
= (aaVc^p^)x(l-l/4.605xLOD) 
where LOD corresponds to the LOD value of the identified QTL. 

Results 

Palm oil fatty acid composition in the £ guineensis cross 
Li\/12Tx DAIOD 

The results for palm oil fatty acid proportions found in the 
LM2T X DAIOD cross are given in Table 1 and Table SI. The 
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Table 1. Means, ranges, variances, and coefficients of variation (CVs) for palm oil fatty acid composition and iodine value in the £ 
guineensis intraspecific cross LM2T x DAIOD. 





Traits 


Mean (n = 88) 


Range 


Variance 


CV ' (%) 


LM2T self Mean (n = 5) 


DAIOD self Mean (n = 5) 


CI 4:0 


0.5 


0.3-1.0 


=0.0 


32.0 


0.2 


1.3 


CI 6:0 


40.5 


32.5-50.0 


9.2 


7.5 


32.9 


45.9 


C16:l 


0.1 


0.1-0.2 


=0.0 


30.2 


0.1 


0.1 


CI 8:0 


5.2 


3.7-8.4 


0.8 


17.1 


6.5 


4.1 


C18:l 


43.5 


35.3-50.0 


7.6 


6.3 


49.9 


35.4 


CI 8:2 


9.3 


5.0-12.2 


1.2 


11.5 


9.0 


12.3 


CI 8:3 


0.2 


0.1-0.4 


=0.0 


21.7 


0.3 


0.3 


C20:0 


0.3 


0.2-0.5 


=0.0 


19.3 


0.3 


0.3 


C20:l 


0.1 


0.1-0.3 


=0.0 


37.0 


0.1 


0.1 


Saturated 


46.6 


41 .0-54.4 


6.4 


5.4 


39.9 


51.7 


Monounsaturated 


43.6 


27.2-50.0 


10.6 


7.5 


50.1 


35.6 


Polyunsaturated 


9.8 


7.1-26.0 


4.0 


20.4 


9.3 


12.6 


Iodine value 


55.3 


49.4-61.2 


4.0 


3.6 


59.4 


52.7 



^CV: Coefficient of variation. 
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principal fatty acids were palmitic acid (C16:0, mean 40.5%) and 
oleic acid (C18:l, mean 43.5%), foUowed by linoleic acid (C18:2, 
mean 9.3%) and stearic acid (C18:0, mean 5.2%). The higher 
coefficients of variation were for the fatty acids present in trace 
amounts (means <1%), such as C14:0, C16:l, C18:3, C20:0, and 
C20:l. The ratio of saturated (46.6%) to unsaturated (53.4%) fatty 
acids was in accordance with the 1:1 ratio, as demonstrated by a 
test (data not shown). 

The mean fatty acid contents of the progeny were equal to the 
mean values estimated for their parents LM2T and DAIOD, 
except for C18:2, as demonstrated by a test (data not shown). 

The normality test (data not shown) showed a normal 
distribution for CI 6:0, CI 8:1, and IV but not for the other fatty 
acids. The histograms (Figure S2) showed discontinuous variations 
for C16:l, C18:3, C20:0, and C20:l, and for this reason, they 
were not considered for further analysis of QTLs. 

In the PC A for the LM2T x DAIOD cross with 16 elementary 
variables associated with agronomic traits and palm oil fatty acid 
composition (data not shown), the first four components explained 
6 1 % of the global phenotypic variation, indicating the existence of 
four groups of correlated traits. The variation explained by each 
component was PCI 21.9%, PC2 38.6%, PC3 51.3%, and PC4 
61.0%. The PCI had loadings mainly for fatty acid traits (C14:0, 
CI 6:0, CI 8:0, and CI 8:1), whereas the PC2 was for traits 
associated with oil palm production (Bwt3_5, Bn6_9, and Bwt6_9, 
Spikelets), the PC3 was for fruit pulp-related traits (%PF, %POP, 
and %KF), and the PC4 was for fruit production traits (Fn, Fwt, 
and %FB). 

The Pearson's correlation coefficients (Table 2) for the most 
relevant fatty acids showed that C14:0 was positively correlated 
with C16:0 and C18:3 and negatively correlated with C18:0, 
C18:l, and IV. Palmitic acid C16:0 was negatively correlated with 
C18:0, C18:l, and IV. Stearic acid C18:0 was positively correlated 
with C18:l and IV. Oleic acid C18:l was negatively correlated 
with C18:2 and C18:3 and positively correlated with IV. 
Considering all vegetative, production, and fatty acid composition 
traits, few correlations were significant at p<0.05 between a fatty 
acid proportion and a vegetative or production trait (Table S2). 
Globally, fatty acid proportions were not correlated with the 



vegetative or production traits under study, i.e. they were 
statistically independent, except for the bunch number and the 
number of leaflets per mature leaf of rank 17. 

Relationships between the five main fatty acid 
proportions in E guineensis and £ oleifera 

The results of the two PCAs for only palm oil fatty acid 
composition in LM2T x DAIOD and SA569 showed that the first 
three components explained 93.2% and 87.3% of the global 
variation in LM2T x DAIOD and SA569, respectively, represent- 
ing high percentages of same order of magnitude (Table 3). The 
factor loadings of the hnoleic acid C18:2 showed that the latter was 
independent of the other main fatty acid traits in the intraspecific 
or interspecific genetic material. However, different tendencies 
were fovmd in the correlations between the five main fatty acid 
proportions, depending on the genetic material. Thus, in the E. 
guineensis cross, the fatty acid proportions of CI 4:0, CI 6:0, C18:0, 
and CI 8:1 were all highly correlated with each other (as mainly 
represented by PCI), whereas in the interspecific SA569, CI 6:0 
and CI 8:1 (associated with PCI) were correlated with each other 
but were independent of CI 4:0 and CI 8:0 (associated with PC2), 
which, in turn, were correlated with each other. 

The projection of the variables in the PCI -PC 2 plans illustrates 
these principle relationships, as shown in Figure 1. 

QTLs involved in fatty acid composition 

Sixteen QTLs associated with palm oil fatty acid composition 
were evidenced by the Kruskal-WaUis (K-W) analysis, with one to 
three QTLs per fatty acid or iodine value (IV) (Table 4 and 
Figure 2). QTLs only detected by the K-W method and 
considered "putative" were mapped in linkage groups (LGs) 1, 
4, and 13. For these putative QTLs, a peak LOD value was 
observed with the IM and/or MQM methods at the same or a 
nearby location, although not significant. Nine QTLs were 
confirmed by IM and ten QTLs by MQM, at the signiflcant 
genome-wide threshold a of 1% or 5%. These latter QTLs were 
located in five LGs (4, 8, 9, 14, and 15) of Eg_Map. 



PLOS ONE I www.plosone.org 



5 



May 2014 | Volume 9 | Issue 5 | e95412 



Genetics of Palm Oil Fatty Acid Composition 

















s 


-0.86^ 


-0.25^ 










o 


0.55' 


0.74' 


0.22 


d 

1 


-0.11 


-0.14 


-0.02 


0.09 


0.13 


-0.12 


-0.07 


* 




* 










-0.40' 


-0.43^ 


-0.45* 


0.50** 


0.34** 


-0.01 


-0.15 


CO 

m 




m 


0.15 


0.25* 


o 




d 


d 


o 


1 


1 


d 




m 


q 


m 


0.04 


0.39* 












1 


1 




















m 


rsi 


CO 
ro 


o 








d 

1 


d 

1 


d 


d 






















o 




oo 
ro 










d 

1 


d 

1 


d 
























o\ 
m 


ro 












d 


d 


























ro 














d 














o 


o 




o 




rsi 


ro 






lb 


CO 


do 


CO 


CO 


u 


u 


u 


u 


u 


u 


U 



O 5 



ON 



< 



O 
-a 



D- 
t/i 
k_ 
0) 

_g 

T3 

ro 

Q 
o 

< 



a. 



o 

J3 



O 
Q- 

E 
o 
u 



'ro 



c 
(U 

c 
o 

Q- 

£ 
o 
u 

"ro 
D. 

'u 
c 



.Q 



o 
< 



d 



I I 



o 


o 


o 




rN 






6d 


do 


6d 


u 


u 


u 


0 


U 



E 



° S 



-o g 

c o 

™ Qj 

-a c 

OJ o 

C Q. 



o ^ 

v 1^ 

g 5 

1- "U 



PLOS ONE I www.plosone.org 



6 



May 2014 | Volume 9 | Issue 5 | e95412 



Genetics of Palm Oil Fatty Acid Composition 



0,25 



-0,2018;; 



LM2T X DAIOD traits (axes PCI and PC2: 84.5%) 

C18L2 




SA569 traits (axes PCI and PC2 : 70.0%) 




18:1 



Figure 1. Principal component analysis (PCA) of palm oil fatty acid proportions of C14:0, C16:0, C18:0, C18:1 and C18:2 in the 
intraspecific cross LiVI2T x DAIOD and in the interspecific pseudo-backcross SA569. Note: the figure show projections on the two first 
axes of the PCA. 
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The percentage of the phenotypic variation explained by a 
significant QTL corrected for the sampling error was low to 
medium and ranged between 10% and 36% (Table 4). The total 
phenotypic variation explained by the QTLs for the principal fatty 
acids was 29% for CI 6:0 (one QTL), 68% for CI 8:0 (two QTLs), 
50% for C 18:1 (three QTLs), and 25% for C 18:2 (one QTL). The 
explained cumulative variation for IV was 40% (two QTLs). We 
did not estimate the LM2T or DAIOD parent effect at QTLs, as 
such an estimation would require a larger mapping population of 
at least 200 palms. 

Several QTLs were closely linked or co-localized. LG 9 showed 
a co-localization of three QTLs for the traits CI 6:0, C18:0, and 
C18:l at the position 0.0 cM and a neighborhood QTL for CI 4:0 
at 9.8 cM. Putative QTLs in LG 4 for C18:0 and C18:2 were 
closely linked along with a QTL for CI 8:1 determined by IM and 
MQM. LG 1 contained closely linked putative QTLs for CI 4:0 
and CI 6:0 and another QTL for IV. Other QTLs for CI 8:1 and 
IV were closely linked in LG 15. Finally, LG 14 and LG 13 
presented one QTL for C18:0 and LG 8 one QTL for C18:2. 

In general, no position correspondence was found between 
intraspecific QTL regions identified on Eg_Map and interspecific 
QTL regions identified previously by Montoya et al. [27] in the 
pseudo-backcross SA569. The two genetic materials showed only 
two cases with QTLs in the same regions: at the bottom of LG 4 
(around mEgCIR0801) and at the top of LG 15, showing C18:l 
and rV QTLs in Eg_Map and CI 6:0 and IV putative QTLs in 
SA569. 

There were two main independent QTL regions, one in the E. 
guineensis cross (in LG 9) and the other in the E. oleifera-denved 
cross (in LG 6), with major effects harboring QTLs for several 
traits with high effects. 



Discussion 

Palm oil fatty acid composition in LM2T x DAIOD 

Our phenotypic data for fatty acid composition show wide 
variation within the cross, with low individual values for palmitic 
acid C16:0 (32.5%) and high values for oleic acid C18:l (50%). 
The mean 1 : 1 ratio of saturated to unsaturated fatty acids was in 
accordance with Ebong et al. [14]). Our data are consistent with 
other reports for genetic materials descending from the La Me 
origin, characterized by its relatively low amount of palmitic acid 
and high amount of oleic acid in comparison to other E. guineensis 
origins [44]. La Me x Deli crosses have concentrations of 40% and 
41% for palmitic acid and oleic acid, respectively [11,16]. More 
recently, Monde et al. [45] evaluated La Me and Deli collections 
from Cote d'lvoire. These authors recorded concentrations of 
3 1 % for palmitic acid and 50% for oleic acid from La Me, whereas 
from Deli, they recorded concentrations of 45% and 38% for 
palmitic acid and oleic acid, respectively. 

The above findings imply that the breeding populations La Me 
and Deli are two sources of variability in fatty acid composition 
and consequentiy provide segregating progenies for evaluating 
genetic variation and searching for QTLs related to the fatty acid 
composition of £. guineensis mature fruits. At the same time, the La 
Me and Deli origins, which are among the major populations used 
by oil palm breeders, represent an important genome resource for 
identifying allelic variants within the E. guineensis species of genes 
involved in palm oil biosynthesis and in determining the final 
proportion of fatty acids in mature fruits. Our LM2T x DAIOD is 
a valuable reference and starting point in that respect. 

[7] and [25] have reported that most fatty acid proportions and 
total unsaturated fats in Elaeis interspecific hybrids are interme- 
diate between the parents' proportions, indicating these are 
quantitative traits with additive effects, with the exception of 
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Figure 2. Sixteen QTLs of palm oil fatty acid proportions and iodine value identified in the £ guineensis cross LM2T x DA10D, 
located on the consensus linkage map in oil palm (Eg_Map) of Billotte et al. [36] and compared to the QTL map for same traits 
published by Montoya et al. [27] in the interspecific £/ae/5 pseudo-backcross SA569. Note: Each microsatellite linkage map has 16 linkage 
groups corresponding to the 16 homologous pairs of chromosomes of the Eloeis genome. The E guineensis Eg_Map {253 loci) is sharing 156 marker 
loci in common and good co-linearity with the linkage map of the pseudo-backcross SA569 (362 loci). The QTLs were identified by the Kruskal-Wallis, 
IM and MQM methods. One star {*) or two stars {**): QTL detected by the MQM method at the genome-wide a threshold value of 5% or 1% 
respectively. No star: putative QTL as only detected by the Kruskall-Wallis test at p<0.005. The names and the positions {cM) of the markers are given 
on the right side of the linkage groups. mEgCIRxxxx and mEgESTxxxx: E guineensis SSR loci. sEgOPGPxxxx: E guineensis gene SNP loci. mCnCIRxxxx: 
Cocos nucifero SSR loci. Marker loci common to both maps are indicated by an extension "_R". The names, positions and confidence regions of the 
QTLs are given on the left side of the linkage groups. In red: are figured the QTLs of saturated fatty acid proportion; in blue: the QTLs of unsaturated 
fatty acid proportion and of iodine value. 
doi:10.1371/journal.pone.0095412.g002 



linoleic acid (CI 8:2), for which E. guineensis seemed to be dominant 
for its corresponding allelic genetic factors. Our study clearly 
confirmed the additive genetic determinants of fatty acid 
proportions in our E. guineensh cross, as the mean values of 
LM2T X DAIOD progeny were the mean values of their parents, 



except for C18:2, which showed a non-additive genetic determin- 
ism, with LM2T appearing dominant over DAIOD, consistent 
with Hardon [7]. 

Similarly, our results on the interspecific pseudo-backcross are 
consistent with Tan et al. [23], who reported a co-dominant 
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heredity in hybrid progenies, with the same exception of linolenic 
acid (CI 8:2), whose effects seem to be dominant in E. guinemsis. 
Our PCA based on the fatty acid variables showed a relative 
independence of linoleic acid C 18:2 compared to other fatty acids. 
This independence might correspond to the de novo fatty acid 
synthesis that occurs in the plastid (C14:0, CI 6:0, CI 8:0, and 
C18:l), under the control of genes with additive effects, whereas 
CI 8:2, a fatty acid whose elongation and desaturation occur in the 
endoplasmic reticulum, might be under the main control of a 
gene(s) with dominant effects. 

Our phenotypic data and those of Montoya et al. [27] showed 
intra- and interspecific variability and mainly additive genetic 
determinism for Elaeis palm oil composition. This finding implies 
the strong possibility of using genetic manipulation to improve the 
unsaturated fatty acid proportions of pabn oil based on the genetic 
values of the genitors to be selected. 

Correlations with production traits 

Our estimated correlations between fatty acid proportions in the 
pulp of mature fruits were similar to those of Noiret and Wuidart 
[17] and Wuidart and Gascon [16] in La Me x Deli crosses and to 
those of Noh et al. [12] in the E. guineemis germplasm collected 
from Angola. There were also similarities in the interspecific 
pseudo-backcross SA569 involving E. oleifera [27] and in a pseudo- 
hybrid [26] . This agreement in findings tends to confirm that such 
individual correlations are valid for the whole Elaeis genus. 

An important new finding of this study is a general correlation 
table between vegetative, production and palm oil fatty acid traits, 
which was never published in oil palm. That latter showed few 
correlations between the palm oil fatty acid composition and the 
elementary production traits in the E. guineemis cross, such like 
between the bunch number at the young age with CI 6:0, CI 8:0 
and CI 8:1. We remark as well similar correlations between these 
fatty acid proportions and the number of leaflets per leave at the 
adult age. As far as we know, no individual phenotypic correlations 
have been reported between pakn oil fatty acid proportions and 
vegetative or production traits of oil palm. Correlations exist 
between fatty acid composition and production in other species, 
such as the olive tree [34,46] and Brassica napus [47]. Our results 
indicate that the palm oil fatty acid proportions of mature fruits 
are globally not correlated with elementary vegetative or 
production traits, while few correlations with some production 
traits should not be ignored. This finding suggests that breeding to 
modify the palm oil composition can be performed in E. guineensis 
without affecting important bunch components, such as the 
percentage of pulp on the fruit and the oil content of the pulp, 
which are key parameters in the elaboration of the final palm oil 
yield. Meanwhile, some cautions should be taken regarding the 
bunch number. According to BUlotte et al. [36], who studied 1,182 
palms from 16 different fuU-sib families, the iodine value was also 
positively correlated (p<0.01) with the percentage of palm oil in 
the pulp (%POP). Here, that correlation was not significant in 88 
LM2T X DAIOD individuals, although our samples had a 
narrower genetic base. Interpretation at the level of the Elaeis 
genus must be supported by further experiments using a large 
panel of various E. guineensis, E. oleifera, and interspecific 
populations. 

Genetic information from mapped QTLs and breeding 
perspectives for palm oil composition 

LM2T X DAIOD was considered as a representative of pure 
intraspecific E. guineensis materials and the interspecific pseudo- 
backcross SA569 as an indirect representative of the E. oleifera 
species. In the latter case, most mapped molecular markers were 



segregated from the interspecific parent SA65T, inherited from E. 
oleifera or (by contrast) E. guineensis grandparent genomes but 
monomorphic in the E. guineensu parent P03228D [27]. A main 
cause of this monomorphism is the high homozygosit)' of 
P03228D, a descendant by selfmg of the DeH DA115D genitor, 
itself descending from only 4 ancestral palms of the Deli origin 
[48]. Moreover, Montoya et al. [27] showed that the detected 
QTLs for fatty acid composition in SA569 were only detected in 
the interspecific hybrid parent (SA65T), i.e. from both Elaeis 
grand-parent genomes while no statistical effect of any E. guineensis 
allele was evidenced in the E. guineensis parent P03228D. 

Comparative genome mapping across species or with other 
genomes has been used primarily to demonstrate events of synteny 
or, conversely, propensity for chromosome rearrangement and it 
has provided valuable insights into the evolution of genomes. Our 
comparison of intra- and interspecific Elaeis genetic maps showed a 
high degree of marker locus co-linearity. Consequently, a common 
set of highly transferable SSRs, with known genome positions, is 
available to search common or specific QTL regions in both Elaeis 
genomes. These types of results are frequendy reported, for 
example, in the genus Ruhus of sub-family Rosoideae [49], in the 
genus Vigna [50], and in Eucalyptus species [51]. 

The comparison of genetic maps of cowpea {Vigna unguiculata L. 
Walp) at the subspecies level (Vigna unguiculata ssp. sesquipedalis) or in 
broader species such as Lotus japonica and soybean had revealed 
differences between taxa but highlighted high conservation zones 
and the syntenic relationships between related crop legume species 
or subspecies [52]. A comparative genetic and QTL mapping 
experiment between white oaks [Quercus robur L. and Q. petraea L.) 
showed a significant number of co-locations for QTLs controlling 
the timing of bud burst. The differences between these species are 
based in the influence of environmental factors or the phenotypic 
plasticity inherent in these species [53]. 

The comparison of QTLs throughout the genetic maps of a 
single genus, in this case Elaeis, provides increased information on 
gene contributions to the phenotypic variation in fatty acid 
proportions in various genetic populations. Common genomic 
regions across populations, involving QTLs of interest, facilitate 
positional cloning and marker-assisted selection of agronomic 
genes. For example, this type of strategy' has been used successfully 
between inter- [Sorghum hicolor L. Moench or S. propinquum) or 
intraspecific [S. hicolor) sorghum populations that showed a high 
degree of marker coUinearity and correspondence for QTL regions 
associated with different traits [54] . 

QTL regions common to both Elaeis genomes are in accordance 
with the alignment of QTL maps in other, related species. 
However, these cases represented only two out of 15 QTL regions. 
This situation is similar to the results of Chen et al. [55] in tomato 
species, where, apart from common genomic positions, 75% and 
85% of QTLs were species-specific for fruit weight and total 
soluble solid content, respectively. 

Some QTLs in this analysis were surely missed or falsely 
identified due to the limited size of our mapping populations [56] . 
The QTL power detection in a cross is the same for common or 
specific QTLs at an equivalent QTL polymorphism. The apparent 
higher frequency of specific QTLs can be explained by parent 
homozygosity at QTLs/genes (therefore undetectable) in regions 
for one but not the other cross. Phenotypic variations might have 
been insufficient to identify some common QTLs, for instance in 
the intraspecific cross. 

In fact, QTLs of different genetic types were revealed. On one 
hand, the oil palm cross identified only intraspecific E. guineensis 
QTLs. Through the heterozygous mapping parents, these QTLs 
are responsible for intra-phenotypic variations of the E. guineensis 
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Figure 3. Effects of £ guineensis versus E. oleifera QTL alleles on the palm oil fatty acid composition, estimated by Montoya et al. [27] 
from the interspecific pseudo-backross SA569. Note: the QTL marker loci were used to perform an ANOVA test (type III, post hoc test of Tukey 
at a = 0.05) to estimate the mean effects of the parent QTL marker alleles on the mean of each phenotypic trait. For the hybrid parent SA65T, the 
species origin of the QTL marker alleles were identified, and the allelic effects at the QTL were therefore estimated by contrast of f. oleifera (grand- 
parent SA49D) against £ guineensis (grand-parent LM2466P). 
doi:1 0.1 371/journal.pone.009541 2.g003 



populations La Me and Deli. On the other hand, the interspecific 
pseudo-backcross had two different types of QTLs: 1) purely 
intraspecific E. oleifera QTLs and 2) interspecific QTLs. The first 
type contributes to intra-phenotypic variations of the E. oleifera 
population. It can share (or not) common genomic regions with E. 
guineensis QTLs. The second type is responsible for between-species 
phenotypic differences. It corresponds to homozygous QTLs (fixed 
genes) in each species, with species-specific alleles responsible for 
between-species phenotypic differences. Such QTLs, undetectable 
in pure E. guineensis or E. oleifera mapping populations, likely make 
up the majority of QTLs detected in the backcross. Complemen- 
tary mapping populations of E. oleifera and respective QTL map 
for fatty acid composition would enable us to compare and further 
validate the genomic regions associated with variations within each 
species, and help in exploring further if the genetic architecture of 
pakn oil fatty acid composition is the same or different between the 
cultivated oil palm and its wild relative E. oleifera. 

Two independent QTL regions, one on the E. guineensis cross (on 
LG 9) and the other on the E. o&j/era-derived cross (on LG 6), can 
be considered most important because they each harbor QTLs for 
several traits with strong effects. This pattern could indicate that 
the phenotypic variability of the palm oil fatty acid composition is 
under the control of different regions of the genomes of these two 
species. However, our results are not sufficient to support this 



hypothesis and more work and analyses of complementary 
experiments should be undertaken before generalizing such 
hypothesis to the whole E. guineensis and E. oleifera species. 

Furthermore, the colocalization of QTLs suggested either 
pleiotropic effects or linkage. This finding indicates that in these 
genomic regions, there is either a single segregating locus affecting 
the biosynthesis of several fatty acids pleiotropically or clusters of 
linked QTLs independently affecting the biosynthesis of the 
different fatty acids. Fine-mapping of these QTL regions and the 
analysis of future sequence data from oil palm will help to 
determine whether linkage or pleiotropy is responsible for this 
colocalization. In E. guineensis, the colocalization of major QTLs 
associated with 01 4:0, CI 6:0, CI 8:0, and C18:l was in agreement 
with the phenotypic correlations we observed between those four 
fatty acids. These results are similar to those reported, for example, 
in jatropha [57] and oat [58]. In parallel, the detection of common 
QTL regions in LG 4 and LG 1 5 tends to confirm that these zones 
could be involved in the same genetic determinism for a portion of 
the variations regardless of the Elaeis species. 

The genetic value of oil palm individuals in terms of the fatty 
acid composition of palm oil can be estimated easily (based on 
measurements on the first bunches produced by the individuals). 
Therefore, the main interest in molecular markers associated with 
fatty acid composition does not lie in marker-assisted selection but. 
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rather, in the possibiUty of optimizing crosses between selected 
individuals to accumulate favorable genes. For instance, this 
process has been successfully implemented to increase the content 
of oleic acid in the peanut [59] . From an operational point of view, 
our study can be considered a first step toward achieving this goal 
in the oil palm, which will require more accurate estimates of QTL 
positions and estimates of their eflfects and, possibly, a QTL 
detection study extended to QTLs with smaller effects. QTLs 
associated with different fatty acids (through either linkage or 
pleiotropy) that colocalize in the same genomic regions will be 
difficult to use if their effects are opposite for saturated and 
unsaturated fatty acids. Obviously, a QTL associated with a single 
fatty acid will be easier to use. When QTLs colocalize in the same 
genomic regions for the two species for a given fatty acid, markers 
specific to the favorable species will have to be developed to follow 
the favorable alleles and accumulate them in E. guineensis. 

We chose which QTLs to compare based on trials under 
difiFerent environmental conditions. Therefore, these QTLs might 
be biased by possible environmental effects on the fatty acid 
composition. However, the high heritabflity of fatty acid propor- 
tions and the location of both experiments in similar agro-climatic 
conditions suggest that the QTLs are accurate and comparable 
between sites. A future genotype x environment experiment should 
test this assertion. 

We have no basis on which to formulate hypotheses on the 
genes underlying the E. guineensis QTLs. The few genes of palm oil 
biosynthesis mapped by Montoya et al. [27] in SA569 were outside 
these QTL regions. These QTL genes should be found by other 
methods, preferably based on the search, sequencing, and 
mapping of genes involved in Elaeis palm oil biosynthesis as well 
as gene expression studies of the pulp of developing fruits in 
various Elaeis genetic materials. These experiments wiU be a next 
step of our research. 

Based on the phenot^pic variability existing between E. guineensis 
and E. oleifera in the profile of fatty acids of their respective oils, we 
hypothesized that the genetic architecture of this trait would differ 
between these two species in terms of the numb(;r and position of 
QTLs and in the phenotypic variance they explain. In the studied 
E. guineensis intraspecific cross and interspecific pseudo-backcross, 
dense microsatellite linkage maps, with a high number of common 
and coUinear SSR marker loci, allowed us to identify and compare 
QTLs. The results of the QTL detection undertakc'n in E. guineensis 
were interpreted in comparison with QTLs for the same traits in 
the interspecific pseudo-backcross. Evidence is a diflercnice in 
terms of mmiber and position of QTL between our intraspecific 
and interspecific mapping populations, which might be only due to 
differences in terms of QTL polymorphism between the studied 
genetic materials, not to presence or absence of coding sequences 
underlying QTL regions in the two Elaeis genomes. The whole 
genome sequences in oil palm and E. oleifera recentiy published by 
Singh et al. [6 1] will greatiy help in solving at this level the genetic 
determinism of the palm oil fatty acid composition. 

Having said, we must consider also different E. guineensis versus E. 
oleifera QTL allele effects on the palm oil fatty acid composition 
(Figure 3), from QTL results of Montoya et al. [27] in the 
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